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Abstract 

In a nanoscale molecular junction at finite bias voltage, the intra-molecular distribution of vi- 
brational energy can strongly deviate from the thermal equilibrium distribution and specific vibra- 
tional modes can be selectively excited in a controllable way, regardless of the corresponding mode 
frequency. This is demonstrated for generic models of asymmetric molecular junctions with local- 
ized electronic states, employing a master equation as well as a nonequilibrium Green's function 
approach. It is shown that the applied bias voltage controls the excitation of specific vibrational 
modes coupled to these states, by tuning their electronic population, which influences the efficiency 
of vibrational cooling processes due to energy exchange with the leads. 
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Nonequilibrium processes in nanoscale quantum systems and particularly quantum trans- 
port processes in molecular junctions [1, 2] have received great attention recently. Coupling 
a single molecule to two metal or semi-conductor electrodes at different electrochemical 
potentials drives it into a nonequilibrium steady-state which can deviate significantly from 
thermal equilibrium. The mechanical stability of such a junction and its transport properties 
strongly depend on the distribution of vibrational energy between the different molecular 
modes. Experiments employing scanning tunneling microscopy (STM) [3-5] or molecular 
junctions [6-8] as well as theoretical studies [9-17] indeed demonstrated nonequilibrium vi- 
brational excitation, as well as conformational changes and chemical reactions induced by 
charge transport under finite bias. 

In this communication we demonstrate that the molecular junction scenario offers a new 
route to mode-selective excitations of molecules. Mode-selective excitations were discussed 
extensively in the past in the context of photo-induced chemical reactions [18-20]. The 
original idea was that energy can be directed into specific mode-excitation or bond-cleavage, 
prior to its statistical randomization due to internal vibrational energy redistribution (IVR) 
processes. Since statistical distribution of energy usually prevails in polyatomic molecules 
and in molecules in condensed phase environments, typically the lower frequency modes 
are excited and the weaker bonds are broken. The challenge of overcoming this statistical 
behavior led to control schemes based on laser excitations [21-23]. While the excitation 
of a molecule with a laser pulse results in transient nonequilibrium dynamics followed by 
relaxation to equilibrium, in a biased molecular junction a stationary nonequilibrium state 
is maintained. 

The external bias voltage, as we will show, can be used to drive the molecule into differ- 
ent nonequilibrium internal energy distributions, thus enabling voltage-controlled selective 
excitations of specific vibrational motions, conformational changes or, in principle, bond 
dissociation. While differences in symmetry between different vibrational modes lead to 
well known selection and/or propensity rules [24, 25] for excitations, in this work we focus 
on bias-controlled selective excitations induced by nonequilibrium effects and use a sim- 
ple model to demonstrate how non-thermal excitations of specific molecular modes can be 
independently controlled by the bias voltage. 

To demonstrate and analyze the possibility of mode-selective vibrational excitation, we 
consider here a representative model for vibronic coupling in an asymmetric molecular junc- 
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tion. The model, schematically depicted in Fig. 1, comprises two electronic states on the 
molecular bridge, which are assumed to be localized at different parts of the molecule. This 
localization can be achieved, e.g., employing electron donating or withdrawing functional 
side groups as illustrated in Fig. 1, and results in an asymmetry in the coupling of the elec- 
tronic states to the leads such that the population of the electronic states can be controlled 
by the bias direction (see below). Furthermore, due to the asymmetric charge distribution, 
the coupling of the electronic to the vibrational degrees of freedom will be asymmetric as 
well, thus allowing for selective current-induced vibrational excitation. 

The Hamiltonian describing vibrationally coupled electron transport in this model is 
given by 

H — Hm + -f^Leads + -f^M,Leads, (l) 

where the different parts correspond to the molecular bridge, the leads, and the molecule-lead 
coupling, respectively. The molecular part of the Hamiltonian, Hu — Hq + H-q with 

Wel A^vib -Wvib _|_ c 

m=l i^=l i/=l * m,m'=l 

involves A^ei electronic states on the molecular bridge with energies and the corresponding 
creation and annihilation operators a}^ and am- The vibrational degrees of freedom of the 
molecular bridge ('internal modes') are described by harmonic modes with frequencies fli,, 
corresponding creation and annihilation operators cj, and c^, and are assumed to be linearly 
coupled to the electronic degrees of freedom. To incorporate vibrational relaxation effects, 
each internal vibrational mode is coupled to a thermal bath described by the Hamiltonian 

JVvib ^bath 

Here, o?^^, dp^ denote the creation and annihilation operator for a bath mode with frequency 
ojp^ , and Tjij^p^ is the system-bath coupling strength. All properties of the bath that influence 
the dynamics of the system are described by the spectral densities Ju{oj) = ^^*'=i vt,i3^^{^~ 
uji3^). In the calculations presented below, we have assumed an Ohmic bath model, i.e. 
J,{uj) = ^uje-'^l'^^'^ . 

The electrodes are modeled as non-interacting electron reservoirs with electrochemical 
potential /xl/r- The electronic states in the right and left electrodes are characterized by 
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their energies e^, and the corresponding electron creation and annihilation operators 6^ and 
bk, resulting in the Hamiltonian i^Leads = XlfceL r ^fe^I^fc- Finally, the molecule-lead coupling 
is given by -f^M.Leads = J2keL;m^^,m^T-,k(^Lbk + ZlfceR,m ^R.^^^R-^'^m^fc + h.c.. Here, VL/R,m are 
dimensionless factors that determine the asymmetry in the coupling of the mth molecular 
state to the right or the left electrode, and the molecule-lead coupling strengths are defined 
as ^L/R,fc and correspond to width functions rL/R(e) = 27r XlfegL/R I^L/R,fcP'^(^ ~ ^k)- In 
the calculation presented below the leads are modeled as semi-infinite tight-binding chains, 
where the width-function for the semieUiptic band is rL/R(e) = 1^^/(2171)2 — (e — /xl/r)^ 
with bandwidth 47 and overall molecule-lead coupling strength ^. 

Nonequilibrium transport in molecular junctions is studied employing two complementary 
approaches, where the technical details are summarized in the Supporting Information [26] 
(see also Refs. 13 and 14). Briefly, the first approach is based on propagating in time a re- 
duced density matrix ps(t), which incorporates the electronic states of the molecular bridge, 
the internal vibrational modes and the vibronic coupling in a non-perturbative manner [27]. 
A master equation for ps{t) is obtained employing perturbation theory in the molecule- lead 
and system-bath coupling as well as the Markov approximation [28-30]. Steady-state ob- 
servables are determined at the infinite-time limit, e.g., the internal mode excitation is given 
by (cj,c,y) = lim^^oo trs[ps(i^)ctc,y]. The second approach is based on nonequilibrium Green's 
functions (NEGF). Specifically, we have applied a method originally proposed by Galperin 
et al. [12], which was recently extended to account for multiple vibrational modes and 
multiple electronic states [13, 14]. This approach is complementary to the master equation 
approach, because it is non-Markovian and non perturbative in the molecule lead coupling, 
while vibronic coupling is treated in an approximate non-perturbative way. Thus, the NEGF 
method provides a more comprehensive description of the interactions between the molecule 
and the leads, while the master equation approach invokes no explicit approximations with 
respect to the electron-vibrational coupling. The good agreement of the results presented 
below demonstrates the validity of both approaches for the selected model parameters. The 
results also show that coherences between the molecular levels, which are not included in 
the master equation approach but in the NEGF scheme, are negligible. 

Various possible mechanisms can result in mode-selective vibrational excitation. In the 
following, we discuss mechanisms which are expected to be particularly relevant in nanoscale 
molecular junctions. In all cases we consider the model described above with two elec- 
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tronic states, which are coupled asymmetrically to the two electrodes and selectively to two 
molecular modes, corresponding to Xu,m,m' = Sm,m'Xu,m, with Ai^2 = -^2,1 = in Eq. (2). 
The specific model parameters are detailed in Table I. We note that the mode-selective 
vibrational excitation reported below is also obtained for models with small to moderate 
off-diagonal vibrational coupling terms, Ai,2 = ^2,1 ~ Ai,i/10, but becomes less efficient for 
larger off-diagonal coupling strengths. 

Fig. 2 shows results for the vibrational excitation of the two internal modes as a function 
of the applied bias voltage for a model with a moderate coupling of the internal modes to a 
thermal bath. The results demonstrate that the vibrational excitation of each specific mode 
depends significantly on the bias voltage and can be controlled by its direction. This finding 
can be rationalized considering that the effective electronic-vibrational coupling (and thus 
the vibrational excitation) of a specific mode depends on the residence time of the electrons 
on the molecule (relative to the vibrational period), and thus on the occupation of the 
electronic state to which the mode is coupled. The inset of Fig. 2 shows that the occupation 
of the two electronic states of the asymmetric junction (see Fig. 1) is sensitive to the bias 
direction, so that in each direction the unoccupied state is associated with lower vibrational 
excitation. 

The mechanism of current-induced vibrational excitation and the dependence on the bias 
(direction) can be significantly more complex and intriguing if the cooling by the external 
thermal bath is less effective and other nonequilibrium effects prevail. As an example. Fig. 
3 depicts results for a model with a weak coupling to a thermal bath. The results show that 
the overall vibrational excitation is significantly larger, revealing efficient current-induced 
heating processes schematically depicted in Fig. 4a. The level of vibrational excitation for 
the two modes depends again significantly and in a nontrivial way on the bias voltage. The 
results presented in Fig. 3b for a model where the two vibrational modes have identical fre- 
quencies, Qi = Q2, and the same vibronic coupling strengths, Ai^i = A2,2, reveal furthermore 
that this behavior is not solely caused by different vibrational excitation energies or vibronic 
coupling strengths but rather results from the asymmetry of the molecular junction. It is 
important to note that the relation between the degree of vibrational excitation of a given 
mode and the occupation of the electronic state, to which the mode is coupled, is reversed 
with respect to the results obtained above with a significant coupling to a thermal bath 
[31] (cf. Fig. 2), i.e., vibrational excitation is less significant for the mode coupled to the 
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occupied electronic state. This intriguing finding suggests that other coohng mechanisms 
are active that are selectively associated with the occupied molecular level. 

Indeed, in the absence of efficient cooling by a thermal bath the vibrational excitation 
of the internal modes due to "heating" processes can be limited only by energy loss to 
the electrons in the leads. One of these cooling mechanisms is provided by the transport 
process itself, where upon electron transmission from one lead to another a vibrational 
quantum is absorbed (Fig. 4b), instead of being emitted (Fig. 4a). Notice that when the 
non-equilibrium vibrational distribution involves significant population of excited states, 
cooling by phonon absorption becomes as efficient as heating by phonon emission. However, 
since transport processes depend on the molecular coupling to both of the leads, they are 
not expected to distinguish between the two electronic states or to selectively influence 
the respective mode excitation. In contrast, cooling processes that involve only one lead, 
directly transfer the junction asymmetry to the vibrational system. Two such processes, 
which involve the creation of an electron-hole pair in the left lead [32] by absorbing a 
single vibrational quantum [5], are shown in Fig. 4c and Fig. 4d. When the applied bias is 
sufficiently large, $ ^ 2ei/2, the asymmetry in the coupling of state 1 and state 2 to the left 
lead (sketched by the thickness of the transmission arrows in Fig. 4c and Fig. 4d) implies 
that most of the electron-hole pair creation processes in the left lead involve state 1. As a 
result of these absorption processes, the mode coupled to state 1 looses significantly more 
energy to the left lead than the mode coupled to state 2. Similar electron-hole pair creation 
processes can occur in the right lead as well. However, provided that the two molecular 
states are located sufficiently above the Fermi energy of the electrodes, electron-hole pair 
creation in the right lead requires absorption of multiple vibrational quanta and is thus less 
efficient. The net result is that for positive bias voltages electron-hole pair creation processes 
in the left electrode provide an exclusive source of cooling for the mode that is coupled to 
state 1, which is the occupied electronic state. Inverting the bias direction, cooling at the 
right lead becomes dominant, and again, the mode coupled to the occupied electronic state 
(mode 2 in this case) is more effectively cooled. The bias direction therefore controls the 
efficiency of the prominent vibrational cooling mechanism in the junction, and selectively 
controls the degree of mode excitation. It is noted that this selective vibrational excitation 
depends on the two electronic states being located sufficiently above the Fermi energy of 
the electrodes, and thus can be switched off by shifting them towards the Fermi energy. 
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One therefore has two possibihties to control this mechanism for mode-selective excitation: 
the direction of the bias voltage, which selects a specific mode, and the gate voltage, which 
controls the degree of excitation of the selected mode. 

In summary, we have demonstrated and analyzed mode-selective vibrational excitation 
in single-molecule junctions induced by nonequilibrium electron transport. The results ob- 
tained here for a generic model show that necessary prerequisites for mode-selectivity in 
molecular junctions involve asymmetry of the molecule-lead coupling and/or the charge dis- 
tribution as well as state-dependent vibronic coupling. This can be controlled in a molecular 
junction by choosing appropriate functional and/or anchor groups. The ability to use the 
external bias in order to direct energy selectively into a given molecular vibrational mode 
suggests the possibility of mode-selective chemistry, e.g. the cleavage of a specific, not 
necessarily the weakest, chemical bond induced by nonequilibrium transport processes. 
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Figure Captions 

Fig.l: Sketch of the model molecular junction for positive (/xl > 0) and negative (//l < 0) 
bias voltages. The top figures depict the population of electronic levels. Occupied levels in 
the left (L) and the right (R) lead are represented by the gray-shaded areas. The population 
of the two electronic states of the molecular bridge, which are represented by two black bars, 
is indicated by filled and empty circles attached to the respective bars, respectively. Strong 
(weak) coupling of an electronic state to a lead is indicated by a thick (thin) arched line. 
The bottom figures give a schematic proposition of a scenario for such a molecular junction. 
Electron withdrawing (i?i, i?2, R3) and electron donating (i?4, i?5, Rq) groups attached to two 
phenyl rings should induce localization and asymmetric charge distribution of the electronic 
states, leading to both bias controlled population of the molecular states and bias controlled 
excitation of the vibrational modes. 

Fig. 2: Voltage-dependent steady state vibrational excitation, demonstrating mode- 
selective vibrational excitation in the presence of a moderate coupling to a "cold" bath 
(/cbT = ImeV). The solid, dashed and dot-dashed hues correspond to increasing vibronic 
coupling. The insets show the respective population of the electronic states. The model 
parameters arc summarized in Table I. The black and colored curves are results of NEGF 
and density matrix calculations, respectively. 

Fig. 3 Voltage-dependent mode-selective vibrational excitation in a model with weak cou- 
pling to a bath, where nonequilibrium effects prevail. The black and colored curves are 
results of NEGF and density matrix calculations, respectively. The results have been ob- 
tained for a model with two vibrational modes with different (a) and identical (b) frequencies 
and vibronic coupling strengths, respectively (see Table I for the parameters values). 

Fig. 4: Mechanisms of electronic-vibrational energy exchange in molecular junctions. 
Mechanisms (a) and (b) represent heating and cooling induced by electronic transport be- 
tween the two leads, respectively. Mechanisms (c) and (d), on the other hand, correspond 
to cooling by electron-hole pair creation at the left electrode. The thin arrows in (d) reflect 
the inefficiency of the process due to weak coupling of the unoccupied state to the left lead 
at positive bias voltages. 
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FIG. 2. 
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TABLE I. Model Parameters (Energy values are in eV) 
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